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ABSTRACT 

The data distributed by the National Space Science Data Center on the geophysical 
parameters of precipitable water, sea surface temperature, and surface-level wind speed, 
measured by the Scanning Multichannel Microwave Radiometer (SMMR) on Nimbus-7, 
are evaluated with in situ measurements between January 1980 and October 1983 over 
the tropical oceans. In tracking annual cycles and the 1982-83 El Nino/Southern Oscil- 
lation episode, the radiometer measurements are coherent with sea surface temperatures 
and surface-level wind speeds measured at equatorial buoys and with precipitable water 
derived from radiosonde soundings at tropical island stations. However, there are 
differences between SMMR and in situ measurements. Corrections based on radiosonde 
and ship data were derived, supplementing correction formulae suggested in the data 
handbook. 

This study is the initial evaluation of the data for quantitative description of the 
1982-83 El Nino/Southern Oscillation episode. It paves the way for determination of the 
ocean-atmosphere moisture and latent heat exchanges, a priority of the Tropical Ocean 
and Global Atmosphere (TOGA) Heat Exchange Program. 
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1. INTRODUCTION 

The ocean and the atmosphere exhibit closely coupled variabilities in the tropics. 
The Tropical Ocean and Global Atmosphere (TOGA) program, organized by the Inter- 
governmental Oceanographic Commission and the World Meteorological Organization, 
is a global study of these phenomena. The TOGA Heat Exchange Project (THEP), 
recently established at the Jet Propulsion Laboratory, is an synergistic attempt to deter- 
mine monthly mean surface heat flux in the tropical Pacific (Liu and Niiler, 1985) using 
satellite data. A major component of the ocean- atmosphere heat exchange is the latent 
heat carried by evaporation. At least three parameters measured by spaceborne sensors 
are required to determine the latent heat flux (Liu and Niiler, 1984); they are surface- 
level wind speed (U), sea surface temperature (T), and precipitable water (W). Precipit- 
able water is the total water vapor in an atmospheric column. The Scanning Mul- 
tichannel Microwave Radiometer (SMMR) on Nimbus-7 measured all three parameters 
during the 1982-83 El Nino/Southern Oscillation (ENSO) episode which had far reach- 
ing economical and ecological consequences (see, e.g., Cane, 1983). The determination 
of surface latent heat flux during this episode using Nimbus/SMMR data is one of the 
priorities of THEP. 

Measurements of Nimbus/SMMR during the first year after launch were evaluated 
by Chang et al. (1984) and Gloersen et al. (1984) with optimistic results. During the 
NASA Sea Surface Temperature Workshops, it was suggested that SMMR calibration 
does not compensate adequately for changes in the engineering environment (see, e.g., 
Milman and Wilheit, 1985). This, together with sensor deterioration, may produce sys- 
tematic errors which are not conducive to the study of short-term climate variability. It 
is, therefore, important for us to examine the error characteristics of SMMR data before 
using them in studying the ENSO episode. The ENSO episode, in turn, provides large 
and clear signals to test the sensor’s sensitivity to interannual variation. 

In this study, the Nimbus/SMMR geophysical data, archived and distributed by 
the National Space Science Data Center (NSSDC), are compared with in situ measure- 
ments at the temporal (monthly) and spatial (2° latitude x 2° longitude) resolution 
specified by THEP. The surface level mixing-ratio (Q) derived from SMMR W using 
the relation by Liu (1986) is also evaluated. As discussed in Section 2, different algo- 
rithms were used to process SMMR data at different periods. It is also important to 
remove error due to algorithm changes by comparing the data with in situ measure- 
ments that have consistent temporal characteristics. Reports from ships of opportunity 
provide the only in situ data set with large coverage and a consistent historical archive. 
Our present knowledge of large-scale fields of U, T, and Q were derived from these ship 
data. They are pivotal to our evaluation and form the basis in the determination of 
correction factors. Ship data are concentrated in the northern part of the tropical 
ocean and are scarce in equatorial and southern oceans. Ship data are, therefore, com- 
plemented by two sets of spot measurements: U and T at equatorial moorings; and W 
derived from radiosonde reports at island stations. 

2. SMMR DATA 

SMMR data were arranged according to the "production year" which starts in 
November and ends in October of the subsequent calendar year. Data from Year 2 to 
Year 5, covering a period from November 1979 to October 1983, were obtained from 
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NSSDC in what were indentified as PARM-LO tapes. The contents and processing 
algorithms for Years 2 and 4 were briefly outlined in the data handbook (NASA, 1985), 
but no documentation is available for Year 3 or 5. The three parameters U, T and W 
were extracted with time, earth location, and relevant data flags. 

Different algorithms were applied to T in different production years to correct for 
errors arising from the temporal variation of engineering environment, but the same 
geophysical algorithm was used to derive U from satellite observations. For Year 2 data 
only, a table of W correction factors that depend on the month and on direction of the 
satellite orbit was provided in the handbook to be applied by the user. Different correc- 
tions for other years had been incorporated in the data distributed by NSSDC (P. 
Hwang, personal communication). 

The U and T data provided by NSSDC do not include areas within 600 km of 
land. Daytime measurements of T are also excluded. Additionally, SMMR was turned 
on only during alternate days. The amount of data available in Year 5 is much less 
than available for other years, due to reasons yet to be identified. 

3. IN-SITU DATA 

Surface Marine Data (Tape Deck 1129) archived at the National Climatic Data 
Center (NCDC) were acquired for 1980-83. Three parameters, U, T, and the dew point 
temperature ( rj, were extracted after the application of a quality control procedure (see 
Liu and Niiler, 1985, for details). Comparisons of SMMR and ship data, described 
hereafter, are based on averages of 2° latitude by 2° longitude monthly bins in which 
there are more than ten ship reports of all three parameters. 

The quality of ship data is unknown. For a quantitative error estimate, coincident 
2° latitude by 2° longitude averages of U, T and Q from ship data were compared. For 
the boreal summer (June, July, and August) and winter (December, January and Febru- 
ary) of 1983, all 2° x 2° areas that had more than 20 ship reports in each of the parame- 
ters per month were identified. The data in each 2° x 2° areas were then divided ran- 
domly into two groups and averaged by month. There were 875 pairs of coincident 
data during the winter and 1002 pairs during the summer. The standard deviation of 
the differences of each pair are shown in Table 1; they approximate the noise of ship 
data at the temporal and spatial scales considered. 

As part of the Equatorial Pacific Ocean Climate Studies, moored buoys were 
deployed on the equator, providing almost continuous measurements of U and T from 
July 1981 to March 1983 at 95°W, and from August 1980 to April 1983 at 110°W, sam- 
pling at every 15 min. This set of data was acquired for evaluation. The average 
height of the anemometer was 3.8 m and the average depth of T measurements was 1 m 
(D. Halpern, personal communication). 

Monthly mean soundings from radiosonde reports at 37 ocean stations, from 1980 
to 1983, were extracted from NCDC Tape Deck 9648. These stations and their loca- 
tions are listed in Table 2 and shown in Fig. 1. Monthly mean precipitable water was 
computed by integration using trapezoidal rule after quality selection (see Liu, 1986, for 
details). 
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Table 1. Standard deviation of the 
differences between coincident pairs 
of 2° latitude by 2° longitude monthly 
averages of ship measurements 


Parameter 

Winter 

Summer 

U (m/s) 

1.15 

mm 

T (°C) 

0.55 


Q (g/kg) 

0.49 

mm 


Table 2. Locations of radiosonde stations used in the study 


Name 

Latitude 

Longitude 

Name 

Latitude 

Longitude 

Pacific Ocean 

Atlantic Ocean 

St. Paul 

57.2N 

170.2W 

Lerwick 

60.1N 

1.2W 

Cold Bay 

55.2N 

162.7W 

Stornoway 

58.2N 

6.3W 

Guadalupe 

29.2N 

118.3W 

OWS L 

57. ON 

20. 0W 

Midway 

28.2N 

177.4W 

OWS R 

47.0N 

17.0W 

Lihue 

22.0N 

159.4W 

Bermuda 

32.4N 

64.7W 

Wake 

19.3N 

166.7E 

Gr. Cayman 

19.3N 

81.4W 

Johnston 

16.7N 

169.5W 

San Juan 

18.4N 

66.0W 

Guam 

13.6N 

144.8E 

Barbados 

13. IN 

59.5W 

Yap 

9.5N 

138. IE 

Curacao 

12.2N 

69. 0W 

Truk 

7.5N 

151.9E 

Trinidad 

10.6N 

61.4W 

Koror 

7.3N 

134.5E 

Ilha Trin. 

20.5S 

29.3W 

Majuro 

7. IN 

171.4E 

Gough 

40.4S 

9.9W 

Ponape 

7. ON 

158.2E 




Atuona 

9.8S 

139.0W 




Pago Pago 

14.3S 

170.7W 




Tahiti 

17.8S 

149.6W 

Indian Ocean 


Nandi 

17.8S 

177.5E 




Norfolk 

29.1S 

167.9E 

Seychelles 

4.7S 

55.5E 

Lord Huwe 

31.5S 

159. IE 

Cocos 

12. IS 

55.5E 

Campbell 

52.6S 

169.2E 

Nouvelle 

37.8S 

77.6E 

Macquarie 

54.5S 

159.0E 

Marion 

46.9S 

37.9E 
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4. PRECIPITABLE WATER AND SURFACE-LEVEL HUMIDITY 

Figure 2 shows the average of precipitable water over global oceans measured by 
SMMR. The values for 1981 and 1982 are quite constant, there is a large annual cycle 
in 1980, and the values drop precipitously starting in June 1983. The 1983 drop was 
probably caused by deterioration of the 21 GHz horizontally polarized channel (P. 
Hwang, personal communication). 

Monthly mean precipitable water from radiosondes ( Vf T ) was compared with the 
monthly mean from SMMR (wj in the 2° latitude by 2° longitude areas in which the 
radiosonde stations are located. The time series for four Pacific stations (Midway, 
Majuro, Atuona and Lord Huwe) are shown in Fig. 3 as examples. At Atuona, the data 
does not start until August 1981. It is shown as an example because its location is 
important with respect to the ENSO episode. In the Pacific, SMMR underestimates W 
for all tropical stations; the difference decreases away from the equator until it becomes 
negligible at midlatitudes. In the high latitude North Pacific (St. Paul and Cold Bay), 
SMMR overestimates W; but in the high latitude South Pacific (Macquarie and Camp- 
bell), SMMR appears to underestimate. In the tropical and midlatitude Atlantic, 
SMMR behaves as it does in the Pacific; but SMMR agrees closely with radiosonde data 
at high latitude North Atlantic stations (OWS L, Lerwick and Stornaway). In the trop- 
ical Indian Ocean, SMMR underestimates at Cocos but not at Seychelles. SMMR agrees 
with the radiosondes at midlatitude station Nouvelle, but overestimates at high latitude 
Marion. Therefore, in the tropics and midlatitudes, there is good agreement on the 
zonal dependence of SMMR errors in the three oceans. The only exception is at Sey- 
chelles but the data at this station may be erroneous. Such zonal dependence is more 
clearly defined in Years 3 and 4. Since, in general, W decreases away from the equator, 
the zonal dependence may reflect the lack of sensitivity of the sensor. At high latitudes, 
the error characteristics are different for the different ocean basins. 

Forward stepwise multivariate regressions were used to determine correction formu- 
lae, for nine time segments, in the form of: 

W r = C x + C 2 W* + C 3 W+ + <? 6 (1) 

where M is the absolute value of the latitude in degrees. The value of the coefficients C 

are listed in Table 3; missing values indicate that they do not pass the 5 % significance 
test. The time segments correspond to one for each SMMR production year except for 
the last year. Starting in June 1983, each month has its own regression due to the rapid 
increase in error. There is consistent significant dependence of the error on W n . The 
adjusted values (W nl ), obtained from (1) in the place of W r , agree more closely with W T 
at all tropical and midlatitude stations except Seychelles. The time series of W nl are 
included in Fig. 3. Equation (1) does not remove all errors. In a few subtropical sta- 
tions (e.g., Wake), W nl is consistently lower than W r during summer months when W 

peaks. The time series are characterized by predominant seasonal cycles. The 1982-83 
ENSO is revealed by both SMMR and radiosonde data as a W deficit east of the date- 
line at Majuro and a W surplus west of the dateline at Atuona. 

Surface-level specific humidity Q n was derived from W nl using the relation by Liu 
(1986) and used to compare with the specific humidity Q t derived from the T d from ship 
reports. Figure 4 shows the temporal variation of Q n and Q, averaged over all 2° x 2° 
areas within ± 30° of the equator that have more than 10 ship reports per month. The 
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Fig. 3. Comparison of the temporal variation of precipitable water vapor derived from 
radiosondes, from uncorrected and from corrected SMMR data, at four stations in the Pacific 





Table 3. Coefficients of SMMR precipitable 
water correction equation 


Period 


01/80-10/80 

11/80-10/81 

11/81-10/82 

11/82-05/83 

6/83 

7/83 

8/83 

9/83 

10/83 


UNCORRECTED SMMR 
CORRECTED SMMR 


TIME (MONTH) 


Fig. 4. Comparison of the temporal variation of the mean 
coincident binned surface-level mixing ratio derived from ship 
reports, from uncorrected and corrected SMMR data, over the 
tropical oceans 




















locations of these areas in January and July 1983 are shown in Fig. 5 as examples. In 
Fig- 4, Q n is consistently lower than Q t during summer and early fall. The reason may 
be an underestimation of W as discussed. However, inaccuracy of the Q-W relation 
from Liu (1986) cannot be completely ruled out. 

For further adjustment of SMMR data, a forward stepwise multivariate regression 
in the form of 

Q. = C,Q n + C 2 T n + C 3 (2) 

was performed on four months (January, April, July, and October) of data for each pro- 
duction year. The coefficients are listed in Table 4. The corrected values (<? ne ) obtained 
in place of Q, in (2) are also shown in Fig. 4. Figure 6 compares the standard deviation 
(STD) of Q n -Q, and Q nc - Q,. The correction (2) reduces not only the mean differences, 
but also reduces the STD by about 2 g/kg. The STD is about twice the values shown 
in Table 1 which approximates the noise level of ship data. The STD approaches the 
root-mean-square (RMS) differences as the mean differences become small. From Table 

4. it is obvious that the SMMR errors have consistent dependence on T. 

5. SEA SURFACE TEMPERATURE 

The global mean T from SMMR data is shown in Fig. 7. The 1983 seasonal peak 
appears to be higher than normal but 1983 is an abnormal year. The temporal varia- 
tion of binned T , and T n averaged over tropical oceans are shown in Fig. 8 similar to 
Fig. 4. SMMR underestimates ship values in general. Again, four months of coincident 
data were used to derive a linear regression in the form of 

T.= C,T n + C 3 (3) 

for each production year. The coefficients are listed in Table 5. The corrected values 
(T nc ) obtained from (3) in place of T, are also shown in Fig. 8. The correction removes 
most of the mean differences. The standard deviation is about 1.1°C, which is about 
twice the ship data noise, as shown in Fig. 9. 

In Fig. 10, T n is compared with measurements at equatorial buoys ( T b ). Only buoy 
measurements within the hour of satellite overpass were averaged. Binned averages of 
T n at both 1°N and 1°S are shown. At 95°W, the meridional gradient is large, as 
confirmed by the sea surface temperature fields operationally produced by the Climate 
Analysis Center (not shown). The temporal variation of T n follows T b at 1°S very closely 
until October 1982, at which time it shifted to follow T b at 1°N more closely. At 
110°W, the meridional gradient is not as large. T n agrees closely with T b at 1°S until 
August 1982. After that, T n is about 1° lower than T b . Correction based on ship data, 
in general, increases T n by about 0.5°. It is known that ship measurements, taken at 
boiler intake, are generally higher than buoy measurements (e.g., Tabata, 1978). Again, 
both SMMR and mooring data show ENSO warming starting in July 1982, superim- 
posed on the regular seasonal cycle. 

6. SURFACE-LEVEL WIND SPEED 

The temporal variation of the monthly mean U n over global oceans is shown in Fig. 
11. Year 5 data appears to peak higher than that of other years. The temporal varia- 
tions of the global mean of coincident binned U, and U n are shown in Fig. 12. U n is far 
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Table 4. Coefficients of SMMR 
surface-level mixing ratio 
correction equation 


Table 5. Coefficients of SMMR 
sea surface temperature 
correction equation 


Period 

<?i 

C 2 

<?3 

Year 2 

0.728 

ITKQ 

-3.011 

Year 3 

0.616 

EEU 

-4.853 

Year 4 

0.614 

0.548 

-7.464 

Year 5 

0.663 

0.474 

-6.486 


Period 


c 2 

Year 2 

1.000 

0.573 

Year 3 

1.037 

-0.339 

Year 4 

1.000 

0.103 

Year 5 

1.013 

0.029 


cn 



Fig. 6. The temporal variation of the standard deviation of the 
difference between surface-level mixing ratio derived from 
uncorrected SMMR data and ship data, as compared with the 
corrected SMMR data 



Fig. 7. Same as Fig. 2, except for sea surface temperature 
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SEA SURFACE TEMPERATURE 



TIME (MONTH) 

Fig. 8. Same as Fig. 4, except for sea surface temperature 



Fig. 9. Temporal variation of the standard deviation of the 
difference between sea surface temperature derived from SMMR 
data and ship reports 









WIND SPEED (m/s) WIND SPEED (m/s) 



Fig. 11. Same as Fig. 2, except for surface-level wind speed 



TIAAE (MONTH) 


Fig. 12. Same as Fig. 4, except for surface-level wind speed 
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higher than U , and has less sensitivity; peak to peak differences are less in U n than in U,. 
The agreement is much improved after application of a correction 

^ = 1.710,-7.5 (4) 

suggested in the handbook except that U nc is higher than U, in Year 5. Temporal varia- 
tion of U nc is also shown in in Fig. 12. Mean differences between U nc and U, are negligi- 
ble except for Year 5; their values determined using four months of data per year are 
shown in Table 6. Application of (4) increases data sensitivity, reduces the bias, but 
also increases the STD, as shown in Fig. 13. The STD of U n -U t is somewhat higher 
than ship data noise. 

In Fig. 14, U ne is compared with buoy wind U b at two locations. Buoy measure- 
ments during the hour of satellite overpass were scaler averaged. U n is an average of 
SMMR wind speed within 1° of buoy locations. The time series are coherent, with devi- 
ations on the order of 1 m/s. Some of the differences may be due to the velocity gra- 
dient in the atmospheric surface layer. SMMR tracks the annual cycles rather well and 
the ENSO event shows up as an intensification of the annual lows starting in October 
1982, particularly at 110°W. 

7. DISCUSSION 

This study is intended only to provide an evaluation of the products released by 
NSSDC to general users with regard to applications for the tropical oceans. It has not 
been an attempt to unravel in detail the sensor problems of Nimbus-7. Comparisons 
between binned satellite and ship data and comparisons between buoy measurements 
and coincident satellite measurements during overpasses were performed. For U, the 
correction suggested by the data source appears to work rather well, except for the last 
year. Further improvement is difficult because the coincident satellite and in situ data 
are highly concentrated around the mean, such that a linear correction could not be 
properly determined. Least square regressions, which minimize the variance of the 
ensembles, would give a small gradient and leave large errors at high and low wind 
speeds. Such a correction, while reducing the STD, desensitizes the satellite data and 
introduces large errors in the bias. The errors of Q n have prominent T-dependence; the 
removal of such n cross-talk " was attempted. Removal of other types of "cross-talk" was 
not attempted because the data did not cover a satisfactory range. 

The Seasat experience suggests that W is the parameter that spaceborne microwave 
radiometers can measure best. The algorithm is relatively simple and the accuracy good 
(cf Alishouse, 1983). In this sense, the Nimbus data, as this study demonstrates, have 
more problems than expected. The corrections suggested in the handbook and those 
incorporated in the data are not adequate even for data prior to the sensor change in 
June 1983. Miscalibration of the 21 GHz horizontal channel is not likely to be the only 
problem. The W derived by Prabhakara et al. (1985), without using the radiance meas- 
ured in this channel, also shows temporal and latitudinal dependence, although the 
error characteristics appear to be different from what we found. The sharp increase in 
error starting June 1983, not found in Prabhakara et al., is very disturbing and is prob- 
ably caused by the deterioration of the 21 GHz horizontal channel. The corrections 
presented here serve only as a temporary measure. A more fundamental correction is 
being sought by F. Wentz for THEP and also by the Nimbus project for general users. 
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Table 6. Correction factors 
for SMMR surface-level 
wind speed 





CORRECTED SMMR 







One of the inherent difficulties of comparing spatially averaged satellite measure- 
ments with spot in situ measurements, as frequently discussed, is the spatial variability. 
Another is due to the vertical gradients in the atmospheric and oceanic surface layer, 
quantitative descriptions of which were given by Liu et al. (1979) and others. While T n 
represents the temperature of the upper few millimeters of the ocean, T b was measured 
at 1-m depth, and T, at depths from 5 to 10 m. Only nighttime SMMR data and buoy 
data are used in this study and, at nighttime, the gradient is small. U b in this study 
was measured at 3.8 m while most ship anemometers are located much higher (20 m). 
Most ship winds are estimated from observations of the sea states through the Beaufort 
scale. SMMR observations are closely related to the sea state through surface emissivity 
but the correction to U n was derived from comparison with buoy data. 

The active participation of scientists studying large-scale ocean-atmosphere interac- 
tion and climate variability in the data management of Nimbus/SMMR is important. 
The algorithm used to exclude data in proximity of land eliminated the opportunity to 
examine some oceanic and climatic features. As there are different levels of data degra- 
dation, the users should be given the choice of using the data. For example, there is no 
data on T n and U n in areas of Micronesia and Polynesia, which are occupied only by 
scattered atolls, but data from these areas are important to the study of ENSO. In this 
sense, all data should be included with those of suspicious quality (close to land) flagged. 
Better documentation and uniform reprocessing of the data would greatly benefit gen- 
eral users. At present, the handbook only gives very brief information on some of the 
data distributed by NSSDC. Corrections have been applied by NSSDC to some data 
but not to all data. Making the corrections suggested in the handbook is no easy 
matter. For example, corrections for water vapor are stratified according to ascending 
and descending orbits but there is no indicator in the data specifying which measure- 
ments belong to ascending or descending orbits. 

The Nimbus/SMMR is sensitive to the annual and interannual variations as 
demonstrated in Fig. 3, 10, and 14. By limiting our evaluation to the tropical oceans, 
we might have bypassed some more difficult problems. A survey of the global fields of 
the difference between satellite measurements and ship data in selected months indi- 
cated uneven distribution in the high latitude North Pacific and North Atlantic that are 
not present in tropical oceans. The validity of SMMR measurements in representing 
large scale mean and anomalous fields in the tropical ocean will be addressed in our 
next report. 
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